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| We present a lattice QCD calculation of the characteristics of the p meson decay. The study is 
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carried out on spatially asymmetric boxes using nHYP-smeared clover fermions in the quenched 
approximation. The resonance mass and coupling constant are calculate using the P-wave scat- 
tering phaseshifts, of the isospin 1=1 two-pion system. We use pion masses m n = 418 MeV and 
m K = 312 MeV. In both cases, the p decay is kinematically feasible. We work on lattice sizes 
N z x 24 2 x 48 with lattice spacing a « 0. 1 fin and N z = 24, 30, 34, 48. 
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I. Introduction 



The study of multi-particle systems is an important step in understanding the strong dynamics 
of hadrons. In this report, we focus on the properties of the p meson decay in the / = 1 two- 
pion system. In the last year, several studies of the two-pion system on fully dynamical QCD 
gauge configurations have been done. These studies have taken steps to achieve a more reliable 
and precise calculation by using: larger variational bases, better stochastic estimators, lower pion 
masses, larger boxes and the inclusion of strange dynamics [1, 2, 3, 4]. All these studies used 
Luscher's method [5, 6, 7, 8, 9], and its extensions to moving frames [10]. In the following, we 
present results for the p decay on 4-D asymmetrical tori in the center-of-mass (COM) frame. 

To study the p decay using current techniques, one must extract the two-pion spectrum in 
the energy regime where the p decays. In order to scan this region, some mechanism to dial the 
relative momentum of the pions is required. In the moving frame formalism, the lattice calculation 
are carried out in a frame with non-zero total momentum. The results are then used to determine 
the two-pion spectrum in the COM. In this framework, the ability to tune the relative momentum 
of the pions is restricted because of the finite step size of the total momentum. In certain systems, 
this restriction may not allow the resonance region to be ideally mapped. A possible solution 
was proposed in which one or two of the spatial directions are elongated [11]. In practice, a 
combination of the two methods would most likely lead to the most efficient calculation. In this 
work, we extract the / = 1 two-pion P-wave scattering phaseshifts in the resonance region. The 
relative momentum of the two-pion system is adjusted by elongating one of the spatial directions. 
We use nHYP-smeared [12] clover fermions with Nf = 2. The calculations are carried out in the 
quenched approximation using pion masses m % = 312 MeV and m % = 418 MeV. We use lattice 
sizes N z x 24 2 x 48 with lattice spacing a w 0.1 fm and N z = 24, 30, 34,48. 

The paper is organized as follows. In Section 2, we discuss the details of the lattice calculation. 
The results are discussed in Section 3, and a summary is given in Section 4. 



2. Methodology 

In order to study the p decay, we need to select appropriate lattice sizes. In this work, we 
consider lattice sizes with spatial volume V s = T]24 x 24 2 where we have introduce an elongation 
parameter T] . The elongation of one spatial direction reduces the rotational symmetry to the the 
tetragonal group [13]. We work in the COM frame and fix the allowed angular momentum of 
the low-lying energies by projecting onto the one-dimensional representation. The projection 
restricts the relative momentum of the two-pion to the set 

Prel Prel = (®A f ° r S ° me W * G Z WMl W * ^ °} ' ^ 2A ^ 

and the invariant energy is 



y/i = 2y/ml+p> (2.2) 

The lowest angular momentum I which couples to the channel is I = 1. As a result, the zero 
momentum mode is excluded. To probe the resonance region, we need to select values of 7] which 
will produce a ground state Eq and, if possible, a first excited state E\ which fall in the region where 
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p decays. In order to do this, we make predictions for the spectral behavior as a function of T] by 
enforcing agreement between the effective range formula 

cotg(j)= , ; 1 P ' (2.3) 
Spun (\-ml) 

and the relation between the scattering phaseshifts [11] 

cot5^) = ^-|^oo(l^ 2 ,J?,l) + ^y|^2o(l^ 2 ,r7,l)|. (2.4) 

where contributions from angular momentum I > 3 have been ignored. In Eq. 2.3, m p is the 
resonance mass and g pnn is the coupling constant defined through the effective Lagrangean 

^eff = gpnnY, e abc(Pl ~ Pl)^ (p)7l b (pi)n c (p 2 ) . (2.5) 
abc 

In Eq. [2.4], q = f| with p defined through the relation 



and 2P\ m are generalized zeta functions whose definition can be found in [11]. To generate an 
expectation for the spectrum, we need to fix g p7Ut , m p and m n . The pion mass is easily determined, 
and the coupling constant is taken to be the physical value g p7r7r =6.06(l) . To estimate the resonance 
mass, we generate an ensemble with 7] = 1. In this case, the ground state is near the resonance mass, 
and we therefore make an initial guess m p = , refer to Fig. 1. Using the estimate for the low- 
lying spectrum, we make a conservative guess for the next value 77. Subsequent values are selected 
by refining the values of g p7l7l and m p by fitting the available data with Eq. 2.3. To extract the 




Figure 1: In the left panel, the estimated energy spectrum is displayed for m n = 418 MeV. The blue dots are 
the calculated energies Eq _1 = 0.48(1) and E^~ l = 0.68(2).The blue lines are the estimated behavior of the 
ground, first and second excited states. They are generated using g pKK = 6.06 and am p = 0.48. The value of 
m p is shown with a dashed line. From the figure, it can be seen that there is a relatively small discrepancy 
between Eq _1 and m p . In the right panel, the final results are shown for two-pion spectrum with m n = 418 
MeV. 
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P M P 

C n (t)= 2^7 + 
p * -p 





C 12 (t) = 





p -p 



Figure 2: Quark contractions for C\\ and Cyi 



low-lying energies, we use the variational method [7]. We construct a 2 x 2 matrix of correlation 
functions 

(2.7) 



C(0 ;7 = £(n|^/W^(0)K U=l,2. 



We use operators 



h(t) = nn(p,t) = ^={n (p,t)n + (-p,t) - 7i + (p,t)7t (-p,f)} 

r)=P3(0,0 = 4^Z{ w ( x 5 073w(x,0-^(x,073^(x,0} 
V2 x 



(2.8) 



where ^(p,/) are interpolating operators for the 7T + and 7t~ mesons with momentum p = ^^3. 
To compute the correlation matrix, C\\ and C12 need to be estimated stochastically. C21 can be 
obtained though the relation C21 = C\ 2 , and C22 is evaluated by computing the quark propagator 



•My,* I 

^source 1 tst 



(2.9) 



where A and B run over the 12 spin-color combinations. In Fig. 2, the quark contractions for C\\ 
and C12 are shown. To compute them, we generate noise vectors with the property 



) — 1 spatial © Ispin © 1 color 



(2.10) 



by placing an independently generated random Z4 noise at each lattice point. In Eq. 2.10, the 
angled brackets (...) represent an average over the noise ensemble. Noise vectors are diluted by 
having support only on a single time slice and for a single spin-color combination. Introducing the 
lattice vectors 

Kc D (x,f|p,f 2 ,§0 = ^5cD(x^|y^ 2 )^ p ' y &(y^2) (2.11) 
y 

VAB^AViMVi.h^k) = £.Wx,f|y,fi) [^ Pl ' y 75%^i|z,^)^ P2 ' z &(z^ 2 )] (2.12) 

where % k represents the k th noise vector, the quark contractions for C\\ and C12 can be constructed 
as shown in Table 2. 
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I e ^ y (vAfl(y,f/|-p,r^p,^&)v^(y,r/|-p,r/|0,ri,§k y 

yAZ? 



I ^' p - y (w A5 (y,/ / |-p,^,^)4 5 (y,^|0,^,^; 

yA5 \ 



I e *' y (vAfi(y^/|p 3 ^|-p J ^ J &)v^(y J f/|-p J f/|0,f| 5 ^ 

yA5 \ 



I <? /p - y (w A5 (y,^i-p,^,&)45(y^/l ^n& 

yA5 



I ^ p - y (w A5 (y,/ t / |p,^,T]^)w^(y,/ L / |0,^,77^) 
yA5 



I (vA*(y,f/|p,*i| -p,^,&)^(y^/|o,^,^)75^ 

yA5 \ 

I (vA5(y^/|-p,^|p,^,^)^(y^/|0,^,T7^)7573 

yAB \ 



Table 1: Stochastic estimators for the quark contractions of C\\ and C12. Here <^ and 77 are independently 
generated noise ensembles and (...) represents an average over the noise. 

3. Results 



In Figs. 4 and 5, we show the effective mass plots for eigenvalues X\ {t,tR) and Aa(t,tR) of the 
matrix M{tfy) = C _1 / 2 (^)C(^)C _1//2 (^) with reference time tR = 6. The two-pion spectrum is 
calculated by fitting X\(tfy) and A2 to the curve 



X(t) = c^ /2 cosh [sTsiNt/l-t)] 



(3.1) 



in the temporal region where the largest energy is dominate. We compute the cross-section and 
phaseshifts from the spectral values through Eq. 2.4. The resonance mass m p and coupling constant 
g p7l7l are determined using non-linear regression with weights wt = <Jj~ l where Oi is the variance of 
the i th spectral value. The width of the resonance is computed using 



r = 



SpTlTl 2 



^p/4-m 2 ) 



3/2 



67T 



(3.2) 



We estimate the physical value Y p h by evaluating Eq. 3.2 with the extracted values of g pnM and 
the physical values of m p and m n . The results are listed in Table 3 and displayed in Fig. 3. For 
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m n [MeV] 


m p [MeV] 


gpnn 


r> [MeV] 


T [MeV] 


312(4) 


924(4) 


5.85(35) 


142(16) 


85(8) 


418(2) 


922(1) 


6.25(27) 


161(14) 


16(2) 



Table 2: Fit results for the resonance mass m p , coupling constant g p%% and decay width T. The physcal 
decay width F p h is computed using the calculated value for g pKK and the physical values of m p and m K . 



the two pion masses, we found coupling constants which agree within error and are consistent 
with the physical value g p%% = 6.06(1). The resonance mass showed no appreciable quark mass 
dependence. The only substantial change was a widening of the resonance at the lower pion mass, 
refer to Fig. 3. Despite working in the quenched approximation, our results are similar to studies 
done on fully dynamical QCD gauge configurations. A reasonable comparison can be made by 
looking at the results of [1], since both studies used 0(a) improved wilson fermions with similar 
lattice spacings and comparable values of m K and m p . In both cases, the coupling constant g p%% 
was consistent with its physical value and did not show significant quark masses dependence. The 
resonance mass m p was also found to have a non-substantial quark mass dependence. This suggest 
that the quenched approximation is able to capture qualitatively, the physical dynamics of the p 
decay. However, a proper comparison would require studies which only differ by quenching effects. 
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Figure 3: Cross-section and phaseshift plots for m n = 312 MeV, upper panels, and m n = 418 MeV, lower 
panels. The dashed line marks the calculated value of m p . In yellow, we shade the region m p — T < a^fs < 

m +r. 
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4. Summary 

In this paper, we discussed the calculation of the p decay on spatially asymmetric 4-D tori. 
By varying the degree of asymmetry, we were able to calculate several well separated phaseshifts 
in the region where p decays. The study was carried out in the quenched approximation and 
showed comparable results with similar studies done on fully dynamical QCD configurations. We 
are currently performing the calculation on fully dynamical QCD configurations. This will give a 
more accurate picture of the quenching effects and allow for a comparison with the moving frame 
formalism used in other studies. 
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Figure 4: Effective mass plots for m K = 312 MeV. The elongation parameter 7] is 1.0, 1.25, 1.42 and 2.0 for 
the upper left, upper right, lower left and lower right panels, respectively. 
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Figure 5: Effective mass plots for m n = 418 MeV. The elongation parameter r\ is 1.0, 1.25, 1.42 and 2.0 for 
the upper left, upper right, lower left and lower right panels, respectively. 
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